Spacetime structures of continuous time quantum walks 
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The propagation by continuous time quantum walks (CTQWs) on one-dimensional lattices shows structures 
in the transition probabilities between different sites reminiscent of quantum carpets. For a system with periodic 
boundary conditions, we calculate the transition probabilities for a CTQW by diagonalizing the transfer matrix 
and by a Bloch function ansatz. Remarkably, the results obtained for the Bloch function ansatz can be related to 
results from (discrete) generalized coined quantum walks. Furthermore, we show that here the first revival time 
turns out to be larger than for quantum carpets. 
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Simple theoretical models have always been very useful for 
our understanding of physics. In quantum mechanics, next to 
the harmonic oscillator, the particle in a box provides much 
insight into the quantum world (e.g. 1 1]). Recently, the prob- 
lem of a quantum mechanical particle initially characterized 
by a gaussian wavepacket and moving in an infinite box has 
been reexamined Surprisingly, this simple system 

shows complex but regular spacetime probability structures 
which are now called quantum carpets. 

In solid state physics and quantum information theory, one 
of the most simple systems is associated with a particle mov- 
ing in a regular periodic potential. This can be, for instance, 
either an electron moving through a crystal 1 5, 6] or a qubit on 
an optical lattice or in an optical cavity L^LalM]. F° r tne elec- 
tron moving through a crystal, the band structure and eigen- 
functions are well known. In principle, the same holds for 
the qubit. However, in quantum information theory, the qubit 
on a lattice or, more general, on a graph is used to define the 
quantum analog of a random walk. As classically, there is 
a discrete fioll and a continuous-time full version. Unlike 
in classical physics, these two are not translatable into each 
other. 

Here we focus on continuous-time (quantum) random 
walks. Consider a walk on a graph which is a collection of 
connected nodes. Lattices are very simple graphs where the 
nodes are connected in a very regular manner. To every graph 
there exists a corresponding adjacency or connectivity matrix 
A = (Aij), which is a discrete version of the Laplace opera- 
tor. The non-diagonal elements Aij equal — 1 if nodes i and 
j are connected by a bond and otherwise. The diagonal ele- 
ments An equal the number of bonds which exit from node i, 
i.e., An equals the functionality fi of the node i. 

Classically, a continuous-time random walk (CTRW) is 
governed by the master equation lfl2l IT3I1 



d 
dt 



Pjk(t) = 



(1) 



where Pjk(t) is the conditional probability to find the CTRW 
at time t at node j when starting at node k. The transfer matrix 
of the walk, T = (Tjk), is related to the adjacency matrix by 
T = — 7A, where we assume the transmission rate 7 of all 
bonds to be equal for simplicity. Formally, this approach can 
be generalized to continuous models like the Lorentz gas 1 14]. 



The formal solution of Eq.Q is 

p jk (t) = (j\e Tt \k). 



(2) 



The quantum-mechanical extension of a CTRW is called 
continuous-time quantum walk (CTQW). These are obtained 
by identifying the Hamiltonian of the system with the (classi- 
cal) transfer operator, H = — T iTll H5L Il6ll . Then the basis 
vectors | k) associated with the nodes k of the graph span the 
whole accessible Hilbert space. In this basis the Schrodinger 
equation (SE) reads 



*-|A)= H|fc), 



(3) 



where we have set m = 1 and h = 1. The time evolution of 
a state \k) starting at time to is given by \k(t)) — U(t, trj)|fc), 
where U(i, t ) = cxp(— iH(t — t )) is the quantum mechan- 
ical time evolution operator. Now the transition amplitude 
ctjk(t) from state \k) at time to state \j) at time t reads 



«ifc(*) = 01 



-iHt 



\k). 



Following from Eq.Q the ctjk(t) obey 



1— a 

dt ■ 



ik(t) = ^Hjiaikjt). 



(4) 



(5) 



The main difference between Eq.(|2j and Eq.® is that clas- 
sically J2jPjk(t) = 1, whereas quantum mechanically 

E<K-fc(*)| 2 = 1 holds. 

In principle, for the full solution of Eqs.Q and (|5} all the 
eigenvalues and all the eigenvectors of T = H (or, equiva- 
lently, of A) are needed. Let A n denote the nth eigenvalue of 
A and A the corresponding eigenvalue matrix. Furthermore, 
let Q denote the matrix constructed from the orthonormalized 
eigenvectors of A, so that A = QAQ -1 . Now the classical 
probability is given by 



Pjk{t) = (ilQe-^Q- 1 ^), 
whereas the quantum mechanical transition probability is 

^ k {t) = \a ih {t)\ i = \m^ K Q- 1 \h)\ 2 . 



(6) 



(7) 



The unitary time evolution prevents that ix^k (t) has a defi- 
nite limit for t — > 00. In order to compare the classical long 
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time probability with the quantum mechanical one, one usu- 
ally uses the limiting probability distribution iflTIl 



Xkj 



lim 



1 f T 

T J dtn i k ^- 



(8) 



In the subsequent calculation we restrict ourselves to 
CTQWs on regular one-dimensional (Id) lattices. Then the 
adjacency matrix A takes on a very simple form. For a Id 
lattice with periodic boundary conditions, i.e. a circle, every 
node has exactly two neighbors. Thus, for a lattice of length 
N, with the boundary condition that node N + 1 is equivalent 
to node 1, we have 



A= (Ay) 







* = J 

i=j±l 

else. 



(9) 



For a lattice with reflecting boundary conditions the adjacency 
matrix A is analogous to Eq.(|9}, except that An = Ann = 1 
and Aim = A^i = because the end nodes have only one 
neighbor. Solving the eigenvalue problem for A, which is a 
real and symmetric matrix is a well-known problem, also of 
much interest in polymer physics I18U19I1 . A different ansatz 
describing the dynamics of a quantum particle in Id was given 
by Wojcik and Dorfman who employ a quantum multibaker 
map |20[. 

The structure of H = 7A suggests an analytic treatment. 
For aid lattice with periodic boundary conditions and 7=1 
the Hamiltonian acting on a state | j) is given by 



H|j)=2|i)-|i-l)-|j + l) ) 



(10) 



which is the discrete version of the Laplacian — A = — V 2 . 
Eg. l llOt is the discrete version of the Hamiltonian for a free 
particle moving on a lattice. It is well known in solid state 
physics that the solutions of the SE for a particle moving 
freely in a regular potential are Bloch functions |5, 6]. Thus, 
the time independent SE is given by 



H|$ fl 



E \<5>6 



(11) 



where the eigenstates | <I?6» ) are Bloch states and can be written 
as a linear combination of states | j) localized at nodes j, 



1 N 



iOj 



|j>- 



(12) 



.7 = 1 



The projection on the state \j) than reads = (j\®e) = 

e~ l6 i /\/~N, which is nothing but the Bloch relation + 
Now the energy is obtained from 

Eas.l ll 11 and d!2t as 



En 



2 cos 9. 



(13) 



For small 6 the energy is given by Eg w 9 2 which resembles 
the energy spectrum of a free particle. 



With this ansatz we calculate the transition amplitudes 
dkj(t). The state |j) is localized at node j and may be de- 
scribed by a Wannier function |5, 6], i.e. by inverting Eq. J12t . 



(14) 



Since the states \j) span the whole accessible Hilbert space, 
we have (k\j) = 5kj and therefore via Eq. J12> also 
($e'|$e) = <W 
Then the transition amplitude reads 



N ^ 



e -iE„t e -i8(k-j) _ 



(15) 



The periodic boundary condition for aid lattice of size N 
requires $ e (iV + 1) = $0(1), thus 9 = 2nTt/N with n g 
]0, N}. Now Eq.(I15> is given by 



a jk (t) 



-at 



i2tcos(2n7r/Af) ~i2-rrn(k-j)/N 



N ^ 



(16) 



For small 6, this result is directly related to the results obtained 
for a quantum particle in a box QEIH because then we have 
E n ~ n 2 . 

In the limit N — > 00, Ea.( ll6t translates to 



lim a kj (t) = — 

N— *oo Z7T 



e ~ i0( " k ~^e i2t cos6 
= ^ fe ^e-^J fc _ J (2<), (17) 



— 7V 

- ik-j-i2t 



where Jt(x) is the Bessel function of the first kind |21]. 
The same result has also been obtained with a functional in- 
tegral ansatz 12211 . From Eq. (ll7t we also see that the first 
maxima of the transition probabilities are related to the max- 
ima of the Bessel function, since we have lirriAr^oo TTkj{t) = 
[Jfc_j(2t)] 2 . However, for an infinite lattice there is no inter- 
ference due to either backscattering at reflecting boundaries or 
transmission by periodic boundaries. 

For higher dimensional lattices the calculation is analogous. 
We note that the assumption of periodic boundary conditions 
is strictly valid only in the limit of very large lattice sizes 
where the exact form of the boundary does not matter |5j,|6Jl. 

Very recently it has been found by Wojcik et al, 12311 . that 
the return probability for a Id generalized coined quantum 
walk (GCQW), which is a variant of a discrete quantum walk, 
has the functional form pkk (tr) = [Jo {2t\fD)] 2 , where r and 
D are variables specified in 12311 . which indeed is of the same 
form as the return probability calculated from Eq.(ll7l>. We 
interpret this as an indication that CTQWs and GCQWs, al- 
though not directly translatable into each other, can lead to 
similar results. However, in l23ll the return probability is cal- 
culated for a particle on a very large circle such that interfer- 
ence effects are not seen on the short time scales considered 
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there. By looking ahead at Fig.^ we see that, indeed, on short 
time scales this is also approximately true in our case of the 
CTQW on the finite lattice. Nevertheless, without going into 
further detail at this point, we note this remarkable similarity 
between CTQWs and GCQWs. 

For a CTQW on a Id circular lattice we calculate the quan- 
tum mechanical transition probabilities TTjk(t). Figure [Ha) 
shows the return probability ir^k (t) for a CTQW on a circle 
of 21 nodes first evaluated in a straightforward way by diag- 
onalizing the matrix A numerically, then by using the Bloch 
function ansatz described above. Both results coincide. For 
comparison we also have computed the return probability for 
the infinitely extended system, see Ea.( ll7> . On small time 
scales all the results coincide. At later times waves propa- 
gating on the finite lattice start to interfere; then the results 
diverge and for a finite lattice one observes an increase in the 
probability of being at the starting node. This happens around 
the time t « TV/2. 




FIG. 1: Plot for a CTQW on a circle of length TV = 21 of (a) the 
return probability and (b) the transition probability to go in time t 
from the starting node to the opposite node on the circle. Time is 
given in units of the inverse transmission rate 7 _1 . The results using 
Eq.Q, long dashed line, and Eq. <16t , full circles, are compared to 
the limit TV — » oo, short dashed line. 

In Fig.^b) the probability to go from a starting node to the 
farthest node on the circle, here to go from node 1 to node 11 
(or 12), is plotted. Again the calculations by the eigenvalue 
method and by the Bloch function ansatz are indistinguish- 
able. As before, also the probabilities for the infinite and for 
the finite systems differ. The difference is more pronounced 
because in time t w TV/ 4 counterpropagating waves from the 
starting node interfere at the opposite node. 

The probabilities to go from a starting node to all other 
nodes in time f on a circle of length TV = 21 is plotted in 
Fig. EJleft). (For a CTQW on a circle the starting node is 
arbitrary.) For small times, when there is no interference, the 
waves propagate freely. After a time t w TV/ 4 the waves inter- 
fere but the pattern remains quite regular. The same holds for 
TV = 20, but the structures are more regular, see Fig. fright). 
This is due to the fact that the number of steps to go form one 
node to another is even or odd in both directions for the even- 




FIG. 2: Plot of the probability for a CTQW on a circle of length 
TV = 21 (left) and TV = 20 (right) over time t to go from a starting 
node to all other nodes. See FigQfor units. 



numbered circle, where it is even in one and odd in the other 
direction for the odd-numbered circle. 
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FIG. 3: Limiting probability distribution \kj for a CTQW on a circle 
of length TV = 20 (circles) and TV = 21 (triangles). 

Figure [3] supports this. The limiting distribution Xkj has 
two maxima for TV — 20, one at the starting node 1 and one 
at the opposite node 11, reflecting the higher symmetry of the 
lattice. Whereas there is only one maximum for TV = 21 at 
the starting node 1 . 

At this point it is instructive to look at very small circles of 
TV = 3 and 4 nodes where the analytic results are still handy. 
With Eq.jl6t we find for the transition probabilities for TV — 
3, 



TTkj (t) 



For TV = 4 we have 




^ cos 3 1 
cos 3 1 - 



I cost 
■ cost 



k=j 



COS t 

sin 4 1 

sin 2 t cos 2 t 



k = j 
k = 2j 
else, 



(18) 



(19) 



where wjj (t) and Wj2j (t) are only shifted by a phase factor of 
7r/2 but equal in magnitude. The limiting probability distri- 
butions are for TV = 3, X11 = 5/9 and X12 — X13 = 2/9 and 
for TV = 4, xn = X13 = 3/8 and X12 — X14 = 1/8, and thus 
support the findings for bigger lattices, e.g. Fig. [5] 

The occurrence of the regular structures is reminiscent of 
the so-called quantum carpets 0001- These were found in 
the interference pattern of a quantum particle, initially pre- 
pared as a gaussian wave packet, moving in a Id box. The 
spreading and self-interference due to reflection of the wave 
packet at the walls lead to patterns in the spacetime prob- 
ability distribution. Furthermore, after some time, the so- 
called revival time, the whole initial wavefunction gets recon- 
structed. For a particle in a box, theses quantum revivals are 
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(almost) perfect and the revival time T follows from the en- 
ergy E n = (nnh/L) 2 /2m = n 2 2Trh/T, where L is the width 
of the boxij. 




node k node k 



FIG. 4: Contour plot of the probability for a CTQW on a circle of 
length N = 21 (left) and N = 20 (right) over longer times t than in 
Fig.[2| Dark regions denote high probabilities. See FigQfor units. 

For very long times, Fig.0]shows a contour plot of the prob- 
ability for a CTQW on a circle of length N — 21 (left) and 
N = 20 (right). There is obvious structure in the interference 
pattern. Furthermore, there are areas on this quantum carpet 
where there is a very high probability, visualized by dark re- 
gions, to find the CTQW at its starting point. Thus, quantum 
revivals also occur for the discrete lattice. However, these are 
not perfect. 

The revival time r is given by Q-kjir) = akj(0). Since 
the transition amplitudes are given as a sum over all modes n, 
see Ea. (ll6> . we cannot give a universal revival time which is 



independent of n. Nevertheless, from Ea.( ll6> we get for each 
mode n its revival time 

T ^ l-cos(2n^) ^T [1 + C ° t2( - /jV)] ' (20) 

where r € N (without any loss of generality we set r = 1). 
From Ea. J20> we find that r„ > t ji+ i for n €]0, N/2] and 
t„ < t ii+ i for n £]N/2, N}. For certain values of n, t„ 
will be of order unity, e.g. for n = N/2 we get r„ = ir/2. 
However, for n << N, Eq.(|20} yields r„ = N 2 /2ttti 2 = 
tq /n 2 , which is analogous to the particle in the box and where 
To is a universal revival time. Thus, the revival times r„ have 
large variations in value. To make a sensible statement about 
at least the first revival time, we need to compare it to the 
actual time needed by the CTQW for travelling through the 
lattice. As mentioned earlier, interference effects in the return 
probability 7Tx,i(t) are seen after a time t ps N/2. The first 
revival time has to be larger than this, because there cannot 
be any revival unless the wave reaches its starting node again. 
Our calculations suggest that the first revival time will be of 
order To. From Fig. 0] we see that the first (incomplete) revival 
occurs for N = 20 at t « 70 > 20 2 /2tt and for N = 21 at 
t re 75 > 21 2 /2tt. 

In conclusion we have shown that CTQWs on regular Id 
lattices show regular structures in their spacetime transition 
probabilities. By employing the Bloch function ansatz we 
calculated quantum mechanical transition probabilities (as a 
function of time t) between the different nodes of the lattice. 
These results are practically indistinguishable from the ones 
obtained by diagonalizing the transfer matrix. We note that 
the results obtained via the Bloch function ansatz can be re- 
lated to recent results for GCQWs. The spacetime structures 
are reminiscent of quantum carpets, but have their first revival 
at later times than what is found for quantum carpets. 
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